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I NTRQDUCT I ON 

This final technical report describes the results of 
work performed at this Laboratory under NASA Grant 
No. NGR-009-033. This Grant was one of the earliest Govern- 
ment supported research projects in the area of amorphous semi- 
conductors. Until recently, solid state physics had been 
concerned primarily with crystalline materials - a concern 
which has led to a revolution in electronics technology. 

The study of the fundamental properties of disordered semi- 
conductors indicates that many of these properties responsible 
for the huge technological success of crystalline materials 
are still present in the disordered materials. In addition, 
the disordered systems hold much promise for contributing their 
own unique attributes to the development of semiconductor tech- 
nology. For example, their insensitivity to structural changes, 
wrought by exposure to hostile ambient conditions such as 
temperature fluctuations and radiation fields, immediately 
suggests their application to certain classes of hardened 
devices . “But , “as“in“the case - of crystalline semiconductor T 
technology, the fundamental physical properites, materials 
characterization and basic theory must establish a firm 
foundation before there can be much chance of significant 
technological advance. 

During the two years covered by this report, we have 
investigated the optical properties of carefully prepared 
amorphous silicon and boron. Much has been learned about 
characterizing these materials and measuring their purities. 

The groundwork has been prepared for understanding the behavior 
and properties of amorphous semiconductors and for the fabrica- 
tion of useful devices based upon such an understanding. A 
number of technical reports have appeared in the scientific 
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literature (see Appendix listing 6 references published with 
NASA support) and a paper on the optical properties of amorphous 
silicon and boron films is being prepared. Our work on dis- 
ordered semiconductors has developed a firm foundation under 
NASA sponsorship and further work is being continued under 
alternate sponsorship. 

This report covers three areas of investigation into the 
properties of non-crystalline materials: (1) optical properties 

of elemental amorphous semiconductors, (2) Mossbauer studies 
of disordered systems, and (3) theoretical aspects of dis- 
ordered semiconductors. These will be described briefly 
below. Background material has been provided in the original 
proposal Detailed results of the completed work have been, 

or will soon be, published. The detailed results which have 
been published are attached to this report. 


Optical Properties of Elemental Amorphous Semiconductors 


Optical absorption studies of solids are a classic means 
of determining fundamental properties of semiconductors such 
as band gaps, density of states, and impurity levels. Optical 
studies have, of course, accompanied and aided the development 
of crystalline semiconducting devices. In the investigation 
of amorphous materials, optical studies are playing a similar 
role, but a great deal of controversy surrounds the inter- 
pretation of the results. The sharpness of the band edge, 
for example, is still not entirely resolved. Some investi- 
gators have observed sharp band edges in amorphous germanium 
and silicon, while others have observed tails. Band tailing 
is expected according to the presently accepted theory of 
disordered systems. 
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Optical studies in amorphous materials are somewhat 
more difficult than similar studies on crystalline systems. 

The samples investigated in this report are in the form of 
thin films. There are problems associated with sample purity 
during fabrication, substrate influence, optical interference 
effects, and surface contamination after the sample is removed 
from the vacuum. Exploring the band tails in amorphous systems 
also requires greater sensitivity and lower background than is 
usually employed for normal crystal studies. 

An optical apparatus was constructed to maintain the 
sample purity and provide the sensitivity necessary to explore 
the spectral absorption of amorphous materials. The optical 
facility developed under the NASA Grant includes a double-beam 
spectrometer for recording transmission spectra over a wide 
range of sample environments. The instrument is extremely 
versatile in that almost any combination of energy sources, 
gratings, and detectors may be used to obtain spectra over 
the range from 2000 A in the ultra-violet to 4/x in the near 
infra-red. The novel electronic detection scheme allows 

-accurate -t ran sm-i-t-t anc e-measurement s. _to be made down to 0 . 1% 

of the reference signal, while the SPEX Model 1400-11 double 
grating spectrometer provides excellent spectral accuracy 
and high resolution. 

The complex refractive index for each sample was 
determined by a computer analysis of the transmission versus 
wavelength data for a pair of similarly prepared thin films 
of different thicknesses. A computer program was developed 
which properly takes into account the multiple reflections 
between the thin amorphous film and the air as well as those 

t 

reflections between the substrate of known refractive index 
(usually quartz) and the film on one side and air on the other. 
Transmission data for two films identically prepared except 
for thickness are analyzed to give the index of refraction n 
and the absorption coefficient a for the amorphous film. 
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Preliminary results for silicon (Appendix I) and boron 
films (Appendix II) have been presented at appropriate technical 
symposia. A comprehensive summary of our findings is being 
prepared for publication in a technical journal with copies 
to be sent to NASA as an addendum to this report. Figure 1 
shows the variation with photon energy of the refractive 
index n and the absorption coefficient a. for similarly 
prepared amorphous and crystalline silicon films. These data 
clearly show the presence of an absorptive tail in the amorphous 
films and the absence of such a tail in the crystalline films. 

The process of crystallization of amorphous films is not 
well understood. We have investigated the crystallization 
process in amorphous silicon films by observing optical 
transmission changes in the films as they are heated at various 
temperatures for different lengths of time. A model for 
describing crystallization process was formulated and the 
silicon film data fitted to the model. We believe the result 
to be applicable to other systems as well as silicon. A 
detailed report appears in the literature (see Appendix III). 


In connection with the optical investigations sponsored 
by NASA, we have for some time been concerned with the pre- 
paration and analysis of ultra pure silicon films. The 
technique of sputter-ion source mass analysis of the films 
has been used to characterize the purity of the films and 
as an aid in obtaining purer films. The technique is also 
useful in determining bonding and clustering of atomic 

neighbors in the films. This work has been described in 

(2) 

earlier Progress Memoranda to NASA . A more detailed 
account of this work, which was not supported by the NASA 
Grant, has been published elsewhere v . The ultra pure 
samples were, of course, used in the NASA optical study. 
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Mossbauer Studies of Disordered Systems 

The Mossbauer technique has been established as a useful 
method for investigating the microscopic structure of matter. 

Its application to the study of disordered solids is outlined 
in the proposal 

Under NASA sponsorship, the experimental apparatus was 
acquired and tested. Efforts were made to design a system for 
investigating amorphous tellurium films by Mossbauer 
spectroscopy. The NASA program was terminated before any 
satisfactory amorphous tellurium films could be obtained. 

Recently we have, however, succeeded in preparing amorphous 
tellurium films which give good Mossbauer spectra. There are 
dramatic differences between the spectra of the amorphous and 
crystallized samples. The experimental work and its interpreta- 
tion are continuing at this Laboratory under alternate sponsor- 
ship (AROD) . Preliminary results indicate that, as we anticpated, 
the Mossbauer technique is extremely sensitive to changes 
in the structure of films. 

Theoretical Aspects of Disordered Solids 

Theoretical work, in collaboration with scientists from 
other institutions, has centered on the extension of the single- 
site Coherent Potential Approximation to include the effects of 
scattering from pairs of sites and off-diagonal randomness. Both 
weak and strong scattering cases have been analyzed. The numer- 
ical results for the density of states have been obtained for 
various disordered binary alloys. Detailed results have been 
published and are included in Appendices IV, V and VI. The 
model has also been applied to investigate the effect of 
magnetic disorder on the spin wave spectrum and the curie tem- 
perature of a randomly disordered spin-| Heisenberg ferromagnet. 

In addition, computer programs for analyzing optical trans- 
mission data were developed and modified for interpreting the 
experiments performed on silicon and boron films. 
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Optical Properties of Amorphous Silicon 
Films .* N. A. BLUM, C. FELDMAN and K. MOOR JAN I, 
Applied Physics Lab. The Johns Hopkins U. — . 
Amorphous films { ~ . 5 to 1 . 0 ny thick) oI ~ well 
characterized pure Si were prepared by vacuum 
deposition on fused silica substrates under 
carefully controlled conditions. Sputter-ion 
mass spectrometry has provided information con- 
cerning the purity and composition of the films. 
Optical transmission studies on films of various 
thicknesses has yielded values for the complex 
refractive index over the wavelength range 
0.4 to 2.5 irp. The absorption spectra clearly 
show the presence of an absorptive tail at the 
longer wavelengths relative to identically pre- 
pared films which were subsequently crystallized. 
The results will be discussed in terms of recent 
models affecting the tailing of the density of 
states in amorphous materials. A double beam 
-spectrophotometer— has— been— designed— for these 
experiments and will be described briefly. 

*Work supported in part by NASA Grant No. 

NGR 21-009-033. 
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Mass Spectrom etry, Optical Absorption and Elect rical 
Properties of Amorphous Eoron Films * 

Chai'les Feldman, Kishin Moorjani and Norman Blum 

Applied Physics Laboratory- The Johns Hopkins University 
8621 Georgia Avenue, Silver Spring, Maryland 20910 

I 

ABSTRACT 

Electrical conduction and optical absorption in pure films 
of amorphous boron are described 1 Sputter-ion source mass spectro- 
metry was used to investigate the purity of the films and showed 
that the previous samples, deposited from graphite crucibles, 
contained a large amount of carbon. Various sources of impurities 
during vacuum deposition were identified. The sputter-ion source 
mass spectrometry also gives information about the boron ion- 
clusters ejected from these films. The cluster spectra shows 
peaks at B 5 , suggesting a coordination number of 5 for amorphous 

f i lms . The temperature depe ndenc e_of— elect rica 1 -resist -i-v i-t-y 

identifies activation energies at 0.14 eV, 9.35 eV and 0.65 eV. 

The optical absorption edge in amorphous boron is considerably 

broadened and shifted toward lower energies as compared to the 
i . ■ 

crystalline case. The absorption coefficient a, for values 

4 -1 2 

^10 cm , follows the expression (hF - Eg) /bp and for lower 

values can be fitted to an exponential The values of the 

parameters E & and Eg are determined to be 0.12 eV and 0.62 eV 

respectively. These results are discussed and compared with 

those reported in the literature. 

♦ 

Supported by Naval Ordnance Systems Command 
Contract N00017-72-C4401 , Task A13B 


I. 


Introduction 


Efforts since the previous Boron Conference have been 
concentrated on the improvement and refinement of techniques for 
producing thin films. Refinements have been made in vacuum tech- 
niques for producing the films and in optical techniques for 
measuring their optical properties. A sputter-ion source mass 
spectrometer has been used to analyze the films from the stand- 
point of purity and ion clusters, and this information has been c 

. ,1 

used to improve the purity of the films. The methods for pro- 
ducing purer films and the conditions during vacuum deposition 
are' described in the next section. The purity of these samples 
as determined by sputter-ion source mass analysis is discussed in 
Section III. It is found that the films described in previous 
reports contain large amounts of carbon. Section IV deals with 
the analysis of boron ion clusters in amorphous films as well as 
in the bulk material. The electrical conduction and the optical 
absorption in the purer samples are discussed and compared with 
earlier results in Sections V and VI respectively. Finally, the 
last section contains some concluding remarks on the present 
findings. 

! . ' . 

II. Film Preparation 

Samples were deposited in a 12" water-cooled stainless steel 
vacuum bell jar. The system was evacuated by means of a mechanical 
pump, an oil diffusion pump and a titanium getter pump. Appropriate 
zeolite and liquid nitrogen traps are included in the system. The 
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pre-deposition pressure was less than 5 x 10 torr. Deposition 

—8 —7 

pressures were between 7 x 10~ and 1 x 10 torr. The deposition 

parameters for the samples described here are given in Table I. 

Sample B-44 was deposited in a carbon crucible at a pressure of 
—6 

2 x 10 in a similar manner to that described in the previous 
Conference . An analysis of this film, to be discussed later, 

showed that it contained a large quantity of carbon ^ Sample B-66 
was deposited from a silver-plated copper, water-cooled crucible, 
and Sample B-62 and subsequent films used water-cooled stainless 
steel crucibles with molybdenum liners. A photograph of this 
crucible and arrangement in the system is shown in Fig. 1. The 
large diameter (3.5 cm) of the source allowed one to obtain a 
fast deposition rate. The elimination of the carbon crucible 
brought about a lowering of the pressure during deposition since 
the outgassing of the graphite crucible was difficult to eliminate 
Substrate temperatures between 100°C and 300°C represented the 
maximum temperature recorded by a thermocouple in contact with the 
back of a substrate dur ing the deposition. Ex cept w hen note d, the 
heating of the substrate was due to radiation from the crucible. 

In Samples 62 and 68, effort was made to deposit at high rates in 
order to limit the occlusion of residual gases. The residual 
atmosphere in the vacuum chamber was monitored with a mass 
analyzer. Typical partial pressures of major gases are shown in 
Table II. 

Bulk boron was obtained from United States Mineral Company 
with the quoted purity of 99.99995. Attempts to crystallize the 
boron films by heating the samples in argon following the deposi- 
tion have been unsuccessful. The films tend to crack and flake 


off their substrate as they are crystallizing. Films deposited on 
substrates maintained at 900° C remain, however, continuous and 
smooth. The X-ray diffra.ction pattern of these films still shows 
diffuse rings. 

9 . 

III. Sputter-ion Source Mass Analysis 

The films were analyzed by means of a sputter-ion source 

mass spectrometer . This spectrometer employs a 10 kev argon ion 

beam which bombards the boron film, ejecting boron ions, boron ion 

clusters, ion complexes, and impurity ions from the film. The 

ejected ions are analyzed in a mass spectrometer. A typical 

mass spectrometer trace is shown in Fig. 2. The analysis of boron 

samples is very similar to that recently reported for silicon 
( 2 ) 

films . Impurity analysis was carried out using an energy 
window of 100 eV and 250 eV. This allowed determination of atomic 
impurities without interference from clusters and complexes. 

The impurity analysis of the films is obtained from the 
mass - spectra - by _ means— of - the“foTTo'wi'ng _ '! approxinfate eqiTa'tTohl 

V / Y x+ 
f x+ - 

v y y b + 

where the I's represent ion contents, and Y's the relative ion 
yield. A knowledge of sputtering yields is thus essential for the 
determination of impurity content. Yields based on elements in 
elemental form are given in Table III, along with impurity analysis 
of the samples listed in Table I. If the impurities are assumed 

* 

GCA Ion Microphobe Analytical Mass Spectrometer (IMS 101B) , 

GCA Corp.,. Bedford, Mass. 


to be in the oxidized state, the impurity content would be 10 - 
15 times less. 

The high content of carbon in Sample l?-44 is clearly due to 

the graphite crucible. This is in contras— to statements found in 

the literature that the vapor pressure of carbon is much lower 

than that of boron at the melting point a:nd consequently one does 

(3) 

not expect many carbon impurities in the c:cron film . The 8% 

carbon impurity in Sample 44 is considerably greater than the 

0.2% equilibrium solubility of carbon in cv_lk boron. If the 

sample were crystallized, the carbon would probably reside on the 

(4) 

grain boundaries as a separate C or B^C phrase . Major impurities 
in the films (see B-68) tend to be silicon and iron. As can be 
seen from the tables, silicon may arrive f—'on the source itself. 

The source of iron is more difficult to cea ermine. Analysis of 
a different piece of bulk boron showed an aron contamination in 
small areas on the surface. The high sili.con, oxygen, carbon and 
iron content of B-62A probably represents the SiC >2 showing through 
in pinholes as these impurities do not appear to the same extent 
in samples deposited on tantalum substraT =;■= . A direct comparison 
between series B-68 deposited on SiC^ (non shown here) and on 
tantalum bears this out. Sample B-71 was deposited on a substrate 
held at 900° C. In this sample the silicon and iron as well as 
carbon and aluminum probably comes from ci.nfusion from the sub- 
strate into the film. The fused silica sv..:-strate contains carbon 
and iron impurities as determined by the same mass analysis 
techniques. Care was taken during this hiarh temperature run to 
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avoid heating any stainless steel since it had been shown that 

/n \ 

iron and chromium vapor could contaminate a sample at 900° C ^ . 

The samples do not appear to be contaminated with molybdenum 
or silver from the deposition crucible. 

The carbon and oxygen impurities in all samples are difficult 

to analyze since there appears to be a surface reaction with the 

sample during sputtering in the mass analysis procedure. The 

carbon and oxygen impurity content in silicon films deposited by 

the same technique and at somewhat higher pressures are considerably 

less. McElligott and Roberts have demonstrated that 0 2 and CO 

are chemisorbed strongly on a deposited boron film and this 

chemisorption appears to be greater than on silicon surfaces (5) . 

The amount of oxygen and carbon occluded in the films during 

deposition can not account for the large quantity of C and 0 

observed. Table II shows the calculated values of the ratio of 

gas molecules striking the substrate/cm /sec , (f ), to the boron 

gas 

2 

atoms striking the surface/cm /sec,(Pg). If all these gas 
.molecules— stick-,— ’ which— is— un-l-i-kei-y- 7 — t-he— sum— of-~the — oxygen^bc~af\Tng 
species, assuming C0 2 and O 2 do not decompose, would give 434 ppm 
oxygen impurity in Sample B- 6 S, while the mass analysis shows a 
6000 ppm in Table III. The figures for 0, C and perhaps N are 
thus too high. 

IV. Boron Ion Cluster Analysis 

The mass spectra shown in Fig. 2 show peaks corresponding to 
B 2 ... Bg. The peaks B 3 to Bg have not been observed in the mass 


r - 
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spectra when a thermal source rather than sputter source is used 
The cluster peaks are complicated by the existence of two isotopes 
of almost equal intensity. A cluster of three boron atoms, for 
example, shows mass peaks at 30, 31, 32 and 33. The analysis of 
the clusters was accomplished by converting all peaks to an 
average main cluster 11 R + peak through the prediction of relative 
intensities. The results are shown in Fig. 3, which shows the 
relative distribution of ion clusters ejected from the boron 
target. Note the deviations in the films from the relatively 
smooth distribution at B+, which represents a cluster of 5 boron 
atoms. In the case of previously studied silicon films, the 
Si^ ion deviated from the smooth distribution Using argu- 

ments discussed previously for silicon, and assuming that the 
distribution is related to short range order, one would draw the 
conclusion from these curves that the films contain a large 
fraction of boron with five nearest neighbors. The short range 

order in the amorphous films may thus be related to the 

L7Jl— _ 

G'-rhombohedraT — phase - as~d escribed by Badzian . The cluster 
distribution in the films appears to differ from that of the bulk 

I / 

I 

crystalline phase, however the analysis of crystallized film 
samples must be completed . before more definite conclusions can 
be drawn. 

V. Electrical Measurements 

Electrical measurements were carried out on the samples 
described here in an effort to distinguish differences between 
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these pure samples and the ones described previously k ' . 

Resistance versus 1/T measurements were made in a pure argon 

atmosphere in a tube furnace with gravity contacts onto titanium 

deposited electrodes. The results are given in Fig. 4 for samples 

containing carbon, samples from a molybdenum crucible, and a 

sample deposited on a hot substrate. In each case, the curves 

could be broken into three straight lines representing three 

activation energies. Above 700° K, all of the samples appear to 

show an intrinsic slope of .65 eV. This would lead to a band gap 

of 1.3 eV, which is exactly that reported previously The 

activation energy at around 0.3 eV and 0.15 eV was also similar 

to those reported previously. There is thus little difference 

between the values of activation energies of more pure films and 

the films described earlier. These values are also similar to 

(9) 

those found by various authors in crystalline boron . The 

1 t . ' ' 

source of the activation energies is unknown, however it should 

be noted that the films may still contain a carbon impurity which 

_ ( 4 ) - — 

could— account - for at least one of the energy levels . 

All of the samples exhibit p-type conductivity and have a 

O A 

resistivity of 10 * 3 - 1CP ohm-cm, measured by both a 2-contact and 
a 4-contact technique. ; The value is considerably lower than . 

10 ohm-cm reported previously by us v , and even lower than that 
reported for crystalline boron . However, this low value of 

resistivity is the same as that found in carefully prepared films 
by "intrinsic thermometer method" . The present films were 

• ° M 7 

deposited at 1000 A/min at less than 10 torr in contrast to the 
previous films deposited at 200A/min at 2 x 10 torr. Also, as 
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discussed in Section III, these improved conditions have led to 
more pure films. The enormous decrease in resistivity with increased 
purity could then only be accounted for by the presence of large 
amounts of compensating impurities in the previous films and 
. perhaps even in crystalline boron. 

-5 

The previous samples deposited at 10 torr would readily 
switch from a low to a high conduction state in the manner previ- 
ously described v . The low resitivity of the present samples 
makes switching difficult since it is hard to impose a high field. 

At liquid notrogen temperature, however, some switching does occur. 

r 

The effects of joule heating and the analysis of pre-switching 

(13) 

non-linear currents has been described elsewhere v . 

VI. Optical Measurements 

The spectrophotometer used to measure transmission coefficient T 
is a double beam instrument in which the monochromatic beam from 
a SPEX Model 1400 double grating monochromator is divided by a 
— beam- splitter and arinwed' tb pass through lcfeTnt ical optical 
elements before being recombined and fo’cused onto the detector. 

One leg of the beam is passed through the sample and the other 
through a reference absorber (usually air) before the beams are 
recombined. The sample and reference beams are chopped at two 
different frequencies (175 Hz and 200 hz) by tuning fork light 
choppers. The two frequency components of the detector signal are 
separated by two lock-in amplifiers in a manner such that a voltage 
proportional to the ratio of the sample signal intensity to the 
reference signal intensity is^ available for display on a strip 
chart recorder. The wavelength resolution of the instrument is a 
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function of wavelength and slitwidth. In all cases it is less 
—3 —3 

than 5 x 10 pm and is usually about 10 pm. . The photometric 
precision is T±0.005 over most of the range of the instrument. 

The photometric accuracy in the visible region, for which 
calibrated standards are available, is also T=0.005. 

The optical constants were determined by measuring the 
transmission through two samples, deposited at the same time, which 
were identical except for thickness. The procedure for obtaining 
the optical constants from the transmission data has already been 
described elsewhere The value of the refractive index n at 

X = 2.5 pm was found to be 3.45, which is to be compared with the 
value 3.1 for crystalline boron. This increase of approximately 
10% is similar to the ones reported for amorphous forms of 
silicon germanium^ 16 ) and many other covalently bonded 

semiconductors which have the same short range order in their 
crystalline and amorphous states. 

T her Absorption coefficient a. as a function of energy in the 
range of 0.6-1. 7 eV is shown in Fig. 5 and compared with earlier 
results on amorphous boron For completeness, the data on 

single crystals of B-rhombohedral boron v ’ ' is also included. 

The two sets of samples used in the present work were from 0.13 
to 0.43 pm thick, allowing accurate measurements of transmission 
up to 1.7 eV. The absorption values in amorphous boron (Fig. 5) 
are higher than in the crystalline boron and the absorption edge 
is considerably broadened and shifted towards lower energies. The 
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situation is 
and germanium 


again similar 
(16,20) 


to 


that in amorphous silicon 


(14,15,19) 


Compared to our previous data, the absorption coefficient 

is somewhat higher in the entire energy region; moreover, the 

sharper rise in absorption starting at 1.3 eV (see Fig. 5) is 

not discernible in the present data. The analysis of the data 

4 - 1 

shows that the values of o;^10 cm , in the energy range 0.9 eV 
— 1.65 eV, can be fitted very well to the expression I 

(hi/ - Eg ) 2 

, Of 

hif 

with the "band gap" E g = 0.62 eV. For lower values of o, the 
data can be fitted to an exponential e hU// '^'a. with the characteristic 
energy E a = 0.12 eV. In terms of the recent models introduced 
to discuss the density of electronic states in disordered semi- 
conductors, Eg corresponds to allowed optical transitions between 
band-like states, while E a is an activation energy between states 

in the exponentially decaying tails of the density of states,, 

induced by disorder. Though such models have been successfully 
applied to explain optical absorption in disordered semiconductor 
alloys and even in amorphous germanium and. silicon, with admittedly 
less success, their usefulness is still a subject of lively 
controversy . 


Conclusions 

In the present paper, techniques for the production of 
relatively pure films of amorphous boron have been described. A 
quantitative method for determining impurity content and ion- 


clustering effects in the films has been discussed. Elimination 
of the graphite crucible and the combination of lower deposition 
pressures and higher deposition rates has drastically reduced 
the carbon content of the films. The reduction in impurity 
content has, however, led to an. enormous drop in the value of 

Q 

resistivity (by a factor of 10~) as compared to the earlier 

2 

films, and the values are even lower (by a factor of 10 ) than 
those reported for crystalline boron. It is therefore postulated 
that previous samples and perhaps crystalline boron contain a 
large amount of compensating impurities. This point needs 
further confirmation by deliberate doping of pure films with com- 
pensating impurities and checking the purity of crystalline boron. 
Such an analysis would also throw light on the effect of compen- 
sating impurities in the switching effect. 

It should be noted that the reduced impurity content has 
not led to any significant changes in the thermal activation 
energies. The values of the activation energies, 0.15 eV, 0.35 eV 
and 0.65 eV are almost identical to those reported, previously . 

This would indicate that the present films still contain the 
impurities responsible for the levels associated with these 
activation' energies but the densities of these levels have been 
reduced. The highest activation energy for conduction (0.65 eV) 
found at temperatures exceeding 700° K would lead to a thermal 
gap of 1.3 eV. However, some caution should be exercised in the 
interpretation of this value or the smaller one obtained from 
optical measurements (0.62 eV) as a "band gap". The relationship 



of these values to the concepts borrowed from band theory of 
crystalline solids is not entirely clear. Such large differenc 
between the values of electrical and optical "band gap” are 
common when dealing with amorphous semiconductors and further 
work, particularly on the annealing of the samples, is planned 
to resolve them. 
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TABLE II 


RESIDUAL ATMOSPHERE DURING BORON DEPOSITION 

SAMPLE B— 68 


GAS 

MASS 

APPROXIMATE 
PARTIAL PRESSURE 
(torr) 


ch 4 

16 

9 x 10“ 8 

1600 

h 2 o 

18 

4x10- 9 

70 

co + n 2 

28 

1x10~ 8 

140 

o 2 

32 

3x 10- 10 

4 

C n H m 

43 

3x10~ 10 
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C0 2 
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2 x 10- 8 

220 
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Figure Captions 


Fig. 1 


Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Arrangement in bell jar showing source crucible 
with molybdenum .liner and electron gun. 

Polyatomic mass spectrum of a boron film 
(Sample 68-C) . 

Relative distribution of ion clusters from 
bulk and three different films of boron. 

Resistance as a function of reciprocal temperatur 
for three different films of boron. 

Optical absorption coefficient as a function of 
energy. 

present data on amorphous boron 

— - previous data on amorphous boron 

(Refr. 1) 

*" data on crystalline boron 

(Refr. 17) 

~ data on crystalline boron 
(Refr. 18) 
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THE CRYSTALLIZATION OF AMORPHOUS 
SILICON FILMS* 

NORMAN A. BLUM and CHARLES FELDMAN 

Applied Physics Laboratory, The Johns Hopkins University, 8621 Georgia Avenue, Silver 
Spring, Maryland 20910, U.S.A. 


Received 31 May 1972 

It is well known that the heating of vacuum deposited amorphous silicon 
films above about 700°C produces an irreversible transformation to the 
crystalline state 1 ’ 2 ). Films deposited on substrates near or below room 
temperature may, furthermore, tend to contain voids 3 ). At higher substrate 
temperatures the relative volume of voids diminishes, but the films may 
begin to crystallize. The lowest practical crystallization temperatures and 
times should be used to avoid introduction of impurities. It is therefore im- 
portant to know in some detail how deposition temperature and subsequent 
annealing influence the approach to crystallinity in amorphous silicon films. 
The crystallization process has been followed by observing optical trans- 
mission changes in the films as they are heated at various temperatures. It 
will be shown that the crystallization process is a gradual one which takes 
place at any finite temperature. Heat treatment which is likely to anneal 
away the voids is also likely to make the sample tend toward crystallinity. 

The samples were prepared by electron beam vacuum deposition onto 
pure fused silica substrates at about 2 x 10" 7 torr at rates from 200 to 300 A/ 
/min with the substrate temperature rising from 200°C to 300°C during the 
deposition. Pre-deposition pressure was about 5x 10 -9 torr. The source to 
substrate distance was about 15 cm. For the experiments described here, the 
film thicknesses were approximately 5000 A. Film thicknesses were measured 
by a multiple beam interference technique and revealed a decrease in thick- 
ness, and thus also in the volume, in going from amorphous to crystallized 
samples, of about 8.2%. Samples were heated in a tube furnace flowing with 
pure argon for a fixed time and temperature and then removed from the 
furnace and examined optically. 


* Work supported by NASA under Grant No. NGR 21-009-033 and by Naval Ordnance 
Systems Command, Contract N00017-72-C-4401, Task A13B. 
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For incident photon energies above about2.0eV both amorphous and 
crystallized Si films are heavily absorbing. At 2.6 eV, for example, crystalline 
films exhibit an absorption coefficient a of about 2x 10 4 cm _1 , while “as 
deposited” (amorphous) films show an a at least three orders of magnitude 
greater. Thus, the transmission of the films at 2.6 eV (4800 A) is a sensitive 
indication of the degree of crystallinity or amorphicity. An arbitraty criterion 
for “onset of crystallinity” used here is that 5000 A thick samples have a 
total external transmittance of 5% (corresponds approximately to a = 6 x 10 4 ) 
at a wavelength of 4800 A. Samples meeting this criterion are nearly crystal- 
lized; further heating results in only a slight increase in transmission. The 
results are independent of the exact details of the criterion as long as it 
gives an indication of a phase state somewhere between the two extremes; a 
criterion based upon maximum rate of change of the observed parameter 
with annealing is clearly most useful. In practice, since the sample trans- 
mission was not monitored during the annealing process, the times are 
approximate. The operational procedure was to estimate the time to reach 
the criterion and anneal for that length of time at constant temperature. If 
the criterion was not met (i.e., the transmission at 4800 A was appreciably 
different from 5% at 4800 A), then the anneal was repeated with another 
sample at the same temperature for a shorter or longer time depending on 
whether the sample was on the amorphous or crystalline side of meeting the 
criterion. The procedure was repeated until a satisfactory result was obtained. 

Using the above criterion, the time to reach “onset of crystallinity” t c was 
determined as a function of annealing temperature T for a sample which 
was divided into pieces, each annealed and analyzed separately. This assured 
that all samples started out identical to one another. The plot of log t z versus 
T~ l , shown in fig. 1, gives a straight line, indicating that the simple rate 
expression 

t c = xexp(E 0 lkT) (1) 

is a reasonable approximation relating t c and T. Above and to the left of the 
line is the area of certain crystallinity, below and to the right lies the area of 
amorphicity. Close to the line is the intermediate region where the sample is 
in transition between the two states. 

In eq. (1), t is associated with the characteristic time of a microscopic 
interaction between neighboring atoms, while E 0 is identified with the activa- 
tion energy between the metastable amorphous state and the stable crystal- 
line state. From fig. 1, t^5x 10 _14 sec and £ 0 ^3.1 eV, both figures are 
reasonably consistent with the identification of t as an interaction time 
between atomic neighbors, and E 0 as an activation energy for self diffusion 
in Si 4 ). 
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CRYSTALLIZATION TEMPERATURE (°C) 

727 679 636 596 560 



Fig. 1. Plot of crystallization time versus the reciprocal annealing temperature, showing 
regions where sample remains mostly amorphous and where it is crystallized. 


The plot in fig. 1 applies only to films prepared between 200 °C and 300 °C. 
The detailed behavior probably also depends on the nature of the substrate 
and other preparation parameters. Films deposited on higher temperature 
substrates show evidence of being already partially crystallized. For a film 
deposited at 500 °C we observed points below those shown in fig. 1 , indicating 
that preannealing had taken place. 

Eq. (1) has a form common to many polymorphic transformations which 
do not include nucleation 5 ). The same short range order prevails in both 
amorphous and crystalline form - both structures contain silicon tetrahedra. 
The transformation consists of a displacement of silicon atoms from meta- 
stable sites over a potential barrier ,E 0 to lower energy crystal sites. In the 
film, this produces a polycrystalline structure with small grain size. 

If the volume of the film in the crystalline phase is F c and the volume 
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transforming is proportional to the untransformed volume, then the rate 
of transformation is 

dV, 

d t-W-V')- 

so that 

K = V[1 - exp (- &„()], (2) 

where V is the total volume, and k 0 is the rate constant, 

k 0 = v exp(- E 0 /kT). (3) 

From (2), the fraction corresponding to the crystallinity condition which is 
untransformed at a time t c is 

V - V 

9 = ~ 7~ C = ex P ( ' - k oQ ■ (4) 


The time t c corresponds to the observed optical criterion for crystallinity 
at a temperature T, and thus determines some untransformed fraction g ; 
this constant fraction may be combined with (3) and (4) to give 

t c = rexp(£ 0 /&T), (1) 

where r=v -1 |logg| and has the interpretation mentioned previously. 

A plot such as that shown in fig. 1 is very useful for experimenters wishing 
to anneal films while: (a) preserving most of the amorphous character; or 
(b) crystallizing the sample without unnecessarily risking contamination or 
physical damage by overtreatment. It should be emphasized that such a plot 
applies in detail only to samples prepared under a given set of conditions; 
for other preparation parameters the slope and intercept of the logt versus 
T- 1 line would be different from that shown in fig. 1. The results show that 
the change from a-Si to c-Si is a gradual one which, to a first approximation, 
may be described by eq. (1). At room temperature the amorphous film is 
stable ; usi ng the experimentally derived constants, the time for crystallization 
of an amorphous film at 300 K is ~3 x 10 33 years! 


The authors are pleased to acknowledge the capable technical assistance 
of Messrs. K. Hoggarth and E. Koldewey in carrying out the work reported 
here. 
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Abstract 

The paper develops a self-consistent method for studying disordered 
systems with diagonal as well as off-diagonal disorder. The method has 
general applicability to any disordered system and in the present analysis 
is applied to a monoatomic system whero disorder arises due to a random 
distribution of vacancies. 


§ 1. Introduction 

The equilibrium and the non-equilibrium aspects of disordered solids 
have attracted a great deal of attention recently, experimentally as well as 
theoretically. (See various papers in the proceedings edited by Mott 
1970.) In particular, many different theoretical approaches are available 
for discussing the properties of disordered binary alloys. The principal 
among these is the recently introduced coherent-potential approximation 
(CPA) by Soven (1967). 

The application of the CPA to a disordered binary alloy is based on 
replacing the atomic potential at each lattice site by an undetermined 
coherent potential. The multiple scattering effects from the actual 
potentials are described via a T-matrix, where the scattering potential is 
the difference between the actual potential and the coherent potential. 
In an exact formulation the coherent potential can be determined self- 
consistently from the condition that the configurational average of the 
T-matrix must vanish. In the CPA this condition is replaced by a weaker 
one requiring the configurational average of fhe atomic T-matrix to vanish. 


f Partially supported by Naval Ordnance Systems Command, Contract No. 
NOw-6 2-0604-c, Task A13B, and NASA Grant No. NGR-2 1-009-033. 

I Supported by National Aeronautics and Space Administration Research 
under Grant No. NASA NGR-09-005^072. 
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The success of the CPA is evidenced by a number of recent papers which 
have used the CPA to discuss the static and the dynamic aspects of 
disordered alloys (Velicky, Kirkpatrick and Ehrenreich 1968, Velicky 

1969, Velicky and Levin 1970, Kirkpatrick, Velicky and Ehrenreich 1970, 
Economou, Kirkpatrick, Cohen and Eggarter 1970, Stroud and Ehrenreich 

1970, Soven 1970). The method, however, is applicable to alloys possess- 
ing diagonal disorder only, thus making it useful mainly for alloys com- 
posed of isoelectronic atoms. Recently, attempts have been made to 
include off-diagonal randomness in approximations similar to the CPA 
(Berk 1970, Foo, Amar and Ausloos 1970). 

In the present analysis a more general self-consistent approach is 
developed which is capable of dealing with solids exhibiting diagonal as 
well as off-diagonal disorder. The general ideas followed are similar in 
nature to those in the CPA. The off-diagonal randomness arises due to the 
randomness in the hopping energy of an electron between nearest-neigh- 
bour atoms. This necessitates the introduction of a wave- vector depen- 
dent coherent potential in contrast to the wave-vector independent co- 
herent potential of the single-site CPA. 

The model and the formalism are described in the next section and the 
results are discussed in the last section. 


§ 2. Formalism 

The present formalism concerns itself with the effect of diagonal and 
off-diagonal disorder on the electronic density of states of a one-band 
system. The method has general applicability but we focus pur attention 
on a monoatomic system. The disorder in such a one-component system 
can be simulated by removing a certain fraction of atoms at random, from 
an ordered lattice. The probability of occupation of a given lattice site 
is then given by c(0 < c < 1), where c = 1 represents the completely ordered 
lattice. 

We write the total Hamiltonian for the above system as a sum of 
Hamiltonians describing the pairwise interaction of a given particle with 
its nearest neighbours. 

^" = 2 ^ ( 1 ) 

a 

where the summation is over all nearest-neighbour pairs a and 

*a=**=Wlpf a l+Wmnflm*^+ W m(*l*Oui + a m *a l ). • ( 2 ) 

Here TV (I = §/x, and W lm = Wn^x m = W ml . The disorder is incorporated 
through the variable /x, which is unity for an occupied site and zero other- 
wise. In the above definitions, e denotes the potential energy at each 
lattice site while W is the hopping energy between two nearest neighbours. 
Finally z is the number of nearest neighbours. 
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The equilibrium quantity of interesst is the electronic density of states 
p{E), given by the well-known expression 

p{E) = — - Im Trace G(E ) (3) 

n 


where the functional dependence of the Green’s function on the complex 
energy E is given by 


G(E) = 


1 

E-Jf 


• ( 4 ) 


and the bar denotes the average over all possible configurations of the 
system. 

To discuss the effects of multiple scattering of an electron from systems 
represented by the Hamiltonian (1), we introduce an effective medium 
described by the Hamiltonian, 

S = 2S 0 (a,*o,+a m *a m ) + S 1 (a,*a m +a m *a,), ... (5) 

where £„( = £„) and ^i( —^im) are respectively the diagonal and the off- 
diagonal coherent potentials, as yet undetermined. It is desirable that 
£ 0 and £ x be determined self-consiatently. 

The multiple-scattering effects are now easily described by the T- 
matrix, 

T = (3f-l.)[l + GT] (6) 

where the configurationally averaged Green’s function is given by 

m~. ....... ( 7 ) 

From the relationship 

G=G + GTG, . (8) 

we note that the configurational average of the T-matrix must vanish. 
That is, 

T = 0 (9) 

The eqn. (9) represents the self-consistent condition which determines 
the coherent potentials £ 0 and £ x ; which in turn determine the density 
of states via eqns. (7) and (3). 

In the above analysis the problem has been treated exactly, but to 
proceed any further some approximations are needed. In the CPA, the 
condition represented by eqn. (9) was replaced by a simpler one requiring 
the configurational average of the atomic T-matrix to vanish. Such a 
single-site approximation was adequate to determine the single, k- 
independent, coherent potential introduced there. Below we discuss the 
two-sites approximation appropriate to the system represented by the 
Hamiltonian (1). 
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The full T-matrix (eqn. (6)) can be written in terms of two-sites 


^-matrices by the conventional expression, 

T =2 t «+ 2 t « a b+ 2 l ^bGt y +, (io) 

a az£f}^y 

where the two-sites ^-matrix is given by 

t a = (^ a -L a )\l + GQ (11) 


In eqn. (10), the restricted summations imply that the successive pair 
indices cannot be equal ; that is, in the third term a ^ and /? ^ y but 
a = y is allowed. 

The self-consistent condition (9) is now replaced by a weaker one, 

K = 0 ( 12 ) 

This matrix equation leads to two equations, one for the diagonal matrix 
elements (t mm = 0) and the other for the off-diagonal matrix elements 
(f lm = 0) and they determine the two unknowns £ 0 and 

In Wannier representation, the matrix elements of the f-matrix are 
written as 


hm—( If Im ~ ^l) + ( l^li — S 0 ) [g 0 t hn + 

+ (W lm -L 1 )[g l t lm +g 0 t 

mm ] (13) 

and 

f mm ~ ( Sq) + ( W m j S 1 ) [g^t + f/ 

+ {W mm -'Zo)[gihm + g<!'mm\> ( 14 ) 

where l and m are nearest neighbours. The quantities g 0 { = G mm ) and 
y x ( = G lm ) are obtained from the matrix elements of the Green’s function, 



exp [ ik . (I — w?.)] 
E — zL 0 — y(k)'L 1 ' 


(15) 


In eqn. (15), the structure factor y(k) = £ exp (ik . A), where A is the 
nearest-neighbour vector. A 

Equations (13) and (14) can now be solved for t mm and t lm and the 
requirement that their configurational averages vanish, leads to the 
equations 




2c(l — c) 


+ ^- 2 {(S 1 2 -S 0 2 )^ 0 -S 0 } 


(16) 
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and 

° = 7T {[" ( w - ^ - (l - S oV V - ( W - 2,) j 

+ ^^{[^-20%! + ^},. (17) 

where 

D, = (<7 „ 2 - gr, 2 ) [( W - S ,) 2 - (i - S 0 ) 2 ] 

+ 2(IF-L 1 ) fl r 1 + 2^-S 0 ) fl ri-l, (18 a) 

D 2 = (g 0 2 - Qi 2 ) ^s x 2 + ^o(z~ S 0 ) J - 22,(7! + 0 - 2S 0 )fif 0 - 1 (18 6) 

and 

^3 = (9o“ — ?i 2 )[Si 2 — S 0 2 ] — 22,j7, — 22^,,— 1. . . (18 c) 

The formulation is now complete since eqns. (16) and (17) can, in 
principle, be solved for the two unknowns 2 0 and 2, in terms of e, W and 
the energy E. The density of states in terms of 2 0 and 2, is easily 
obtained by rewriting eqn. (3) as 

p(E)~ — — Im G mm (E) = - — Im <7 0 (19) 

IT 7 T 

or 

Pat.(#)= -- Im?o. ( 2 °) 

7T 

where p at (E) is the density of states/atom and g 0 is obtained from eqn. (15). 

The numerical solution of eqns. (16) and (17) for a three-dimensional 
solid is not entirely simple. The corresponding calculations are now under- 
way and will be reported in the near future. However, the simpler prob- 
lem of one-dimensional disorder chain lends itself to an exact analytic 
solution in the present formulation and is discussed in the next section. 

§ 3. Results 

For the case of one-dimensional chain, eqns. (16) and (17) can be 
decoupled to obtain a quartic equation for p = {2Wg 0 )~ l , written below in 
the form of a dispersion equation, 

1 r c 2 c 2 2c(l — c) 2(1 -c)n 

-- 7 + T + — -+— U = 0, . (21) 

p \_p+x - 1 p + x + 1 p+x p+x+S J 

where x= E — e/2W and 8 = e/2 IF. In these units, the band for an 
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ordered chain (c = 1) is centred at the origin with the half-bandwidth 
equal to unity. 

The numerical values of the density of states/atom (actually 
p at ' = 7r2 W ) obtained from the solution of eqn. (21) are plotted in the 
figure. For ordered lattice (c=l), p(x) has the well-known symmetric 
shape with singularities at the band edges (x = ± 1 ). As soon as c deviates 
from unity the singularities disappear, the main band narrows and a band 
associated with the vacancies appears below the main band. As expected 
(figure), the increasing value of S leads to the shift of the centre of the 
‘ vacancy ’ band away from the lower edge of the main band and also 
causes its broadening. The broadening also results from increasing 
disorder represented by decreasing value of c. 


Fig. 1 




The density of states/atom (in units of 2ttW) versus normalized energy 
x = E — €j2W. The scale on the right refers to the main band and that 
on the left to the ‘ vacancy ’ band. ■ 





Conduction in Low-mobility Materials 


173 


The important question of whether these states, for a three-dimensional 
disordered solid, are localized in nature, can only be answered by calculat- 
ing conductivity. Such a calculation, in the formulation presented here, 
is now under way. 
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COHERENT POTENTIAL THEORY OF 
OFF-DIAGONAL RANDOMMESS: BINARY ALLOY 
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K. MOORJANI** 
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and 

S. M. BOSE* 
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Pennsylvania 19026, U.S. A. 

rhe single-site coherent potential approximation (SS-CPA) for a disordered binary alloy 
s extended in a self-consistent manner to the case of off-diagonal randomness. The density 
jf states for a bcc lattice is calculated in the split band limit for no correlation between 
diagonal and off-diagonal randomness and compared with the ordered and SS-CPA density 
af states. 

1. Introduction 

Ever since Soven 1 ) introduced the single-site coherent potential approxi- 
mation (SS-CPA) for the study of disordered systems, many papers 2 ) have 
applied the approximation to the calculation of the properties of disordered 
solids. The SS-CPA is applicable only to the case of diagonal randomness, 
and these results have been criticized by Stern 3 ). 

Previous attempts have been made to extend the CPA to oflf-diagonal 
randomness 4 ). A self-consistent extension is discussed here. 


2. Theory 

The disordered binary alloy A t Bj_* is characterized by the Hamiltonian 
//=r ( |/) E ,</|+T I ^|/>JF (m <m|, (1) 

* Supported in part by NASA grant # NGR-09-005-072. 

** Supported in part by Naval Ordnance Systems Command, Contract # NOw-62-0604-c, 
Task A13B, and NASA grant #NGR-2t-009-033. 

t Present address: Harry Diamond Laboratories, Washington, D.C. 20438, U.S. A. 
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where the atomic energy, e,, can take the value e A or e B and the off-diagonal 
matrix element, W lm , can assume the values W AA , W AB , or W BB . 

In the spirit of the SS-CPA, a configuration-independent Hamiltonian is 
introduced by 

H 0 =z l \iyz 0 o\ + i,*JJ>r 1 <m|, (2) 

where X 0 = 2' // and I l =Z lm are the diagonal and off-diagonal coherent 
potentials, respectively. The index m in each of the Hamiltonians will be 
restricted to the first nearest neighbor of /. The Green’s functions satisfy the 
equations (E—H) G = 1 and (E— H 0 ) G 0 = 1 and are related by G = G 0 + G 0 x 
x TG 0 , where T = H— H 0 and the scattering T-matrix is given by T =T + EG 0 T. 
The self-consistency criterion, <G> = G 0 yields <T)=0, where the angular 
brackets indicate a configurational average. 

Since the operator T consists of diagonal and off-diagonal components, 
it is convenient to introduce two T-matrix equations corresponding to the 
scattering from the two components, i.e., 



Ti — Tj + T iG 0 T t , i — 1,2, 

(3) 

where 

=z l r(i) = z,\iy(e,-z 0 )(i\, 

(4) 

and 

r 2 = Zi* m r (/, m) = Z,* m \ /> (W lm - 2T 1 )<m| . 

(5) 

The total J-matrix can then be written as 



T = T l + T 2 + (T,G 0 BF + T 2 G 0 AU + AU + BV), 

(6) 

where 

A = T t G 0 T 2 , 

(7) 


B — T 2 G 0 T 1 , 

(8) 


U =(1 — G 0 A )~ l , 

(9) 


V = (l — G 0 B)~ l . 

(10) 

The first two terms in eq. (6) represent independent scattering from r l 

and r 2 . 


while the last terms enclosed by parenthesis represent the correlated scatter- 
ing of r l and r 2 . 

The two operators, r l and T 2 , are written as sums over single entities, i.e., 
over single site and pair operators as in eqs. (4) and (5). Then 

r, =T i r(0 + r,„, 2 r(/ I )G o r(/ 2 ) + - (ii) 

This can be written in terms of the single-site /-matrix 

/(/) = r(/) + r(/) g 0 t (/), 

T t =Z l t(l) + Z h * h t(l i )G 0 t{l 2 ) + -.. 


as 


( 12 ) 
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Since there are two coherent potentials, only two matrix elements of T are 
needed, viz., T„ and T lm , l and m being first nearest neighbor pairs. Simul- 
taneously, intermediate states appearing in these matrix elements are restrict- 
ed to first nearest neighbor pairs. A set of terms consistent with this pair 
approximation must be extracted from the /-/ and l-m matrix elements of 
eq. (12). Hence, scatterings from / and its nearest neighbor site must be con- 
sidered. The diagonal elements of <7)) in the pair approximation are 

<(!))„> = <h> + z {t 2 g \tj{ 1 - glhQ-) , (13) 

where /, = (/|/(/)|/), g 2 = (l\Go\fn) and Z = number of first nearest neighbors. 
Eq. (13) represents the SS-CPA /-matrix vertex corrected for scattering from 
its first nearest neighbor. 

It can be shown that the off-diagonal matrix elements of < T x > are 

= Z <<7ohhn/(l — ■ (14) 

If scattering off the nearest neighbor atom is neglected, i.e., g ,=0, then 
<( 7 i)u> = <^> and <(T 1 ) im > = 0 which is the SS-CPA result. 

Now, if a =(/, m) designates an (/, m) pair, then 

r 2 = z x r (a) + z aif r (a) G 0 r(p) + -, (15) 

which can be written in terms of a pair-site /-matrix, /(a)=T (a)+E (a) G 0 x 
x/(a), as 

T 2 =I a f (a) + r a ^t(a) G 0 (()?) + •••. (16) 

Then in the pair approximation 

<(T 2 )„>=Z</„>, (17) 

and 

<(r 2 ) (m > = z</ lm >, (is) 

where /n = (/|/J/) and /, m = (/|/J m ). The matrix elements, /„ and t lm , are 
found by taking the diagonal and off-diagonal matrix elements of the pair- 
site /-matrix and solving the resultant set of coupled equations. Then 

tn = r 2 ml 9 ol{{\-r migi ) 2 - r 2 ml g 0 ), (19) 

and 

hm — (1 — r, i m g i )/(( 1 — Eml^l) 2 — (20) 

Given the matrix elements of 7) and T 2 are known, the diagonal and 
off-diagonal contributions to the correlation terms in eq. (6) can be found. 
The self-consistency condition, < T„> = < T (m >= 0, results in two non-linear 
equations for two unknowns, Z 0 and 2). 



O'Ol 


158 


T. TANAKA ET AL. 



ENERGY 

Fig. 1. The density of states (in units of tt Z W AA ) versus reduced energy {E-e a )/Z W AA for a bee lattice. Wbb = W-aa, 
fV A B = ytV A A, X = concentration of A atoms and <5 = (6a — sb)/Z W \ a , where £b= — e A . The scale on the right is for the 
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3. Results 

These two equations without correlation were solved iteratively for a bcc 
lattice. The density of states is shown in fig. 1 and compared with the SS- 
CPA. The incompleted part of the curve in the impurity band results from 
the lack of convergence in the solutions. The novel feature is the appearance 
of structure in the impurity band which is absent in the SS-CPA. Further 
details and discussions will be published elsewhere. 
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A one-band model of a random binary alloy A B, 

X X— X 

is analyzed in terms of a two-sites coherent potential approx- 
imation. In the tight-binding Hamiltonian,, the off-diagonal 
randomness is introduced via the composition dependent 
hopping energies between nearest neighbor sites. The in- 
clusion of the off-diagonal randomness correlates the 
scattering from a given site to that from its nearest 
neighbors. Such a correlation is incorporated in the 
handling of diagonal randomness (arising from the composi- 
tion dependence of the atomic potentials) by treating the 
diagonal randomness in the pair approximation. The theory 

leads to the wave-vector-dependent coherent potentials and 

• 

previous approximations used in this problem are easily 
obtained in appropriate limits. The numerical results for 
the electron density of states are presented for a number 
of different alloys and compared with earlier calculations. 

For the case of diagonal randomness only, the present 
theory results in the appearance of structure in the 
density of states of the minority component band. This is 
in contrast, to the results obtained from the single-site 
coherent potential approximation, but in agreement with 
the recent work of Schwartz and Siggia. The presence of 
the off-diagonal randomness leads to further structure in 
the density of states. 



I. 


INTRODUCTION 


The coherent potential approach, based on the 
multiple scattering formalism of Lax (1951), has proved 
to be a useful method for the investigation of disordered 
alloys. A simple and elegant version of this approach, 
referred to as the single-site coherent potential approxima- 
tion (SS-CPA) (Soven 1967, Velicky et al. 1968) has been 
fruitfully applied to a random binary alloy A B in which 

X 1“X 

disorder arises only due to the difference between the atomic 
potentials of the two components of the alloy. The hopping 
integrals in the tight-binding Hamiltonian are assumed 
to be independent of the composition of the alloy and hence 
are translationally invariant. Thus, only the diagonal 
part of the Hamiltonian is assumed to be random. 

The SS-CPA has been elucidated (Velicky et al. 1968) 
and applied to various semiconducting and metallic alloys 
(Stroud and Ehrenreich 1970, Levin and Ehrenreich 1971, 
Economou et al . 1970). It has been . extended to the calcula- 
tion of transport coefficients (Velicky 1969, Levin et al. 
1970) and optical absorption (Velicky and Levin 1970) and 
its equivalence to previous approaches based on atomic 
picture (Matsubara and Toyozawa 1961, Yonezawa and Matsubara 
1966, Matsubara and Kaneyoshi 1966) (in contrast to the 
CPA, which is based on strating from an averaged crystalline 
solid) has been demonstrated (Leath 1970, Ducastelle 1971). 
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However, a recent criticism of the SS-CPA is 
worth noting. As Stern (1971) has pointed out, the numer- 
ical work based on the SS-CPA has very little applicability 
to real alloys. The essential weakness of the model lies 

in the perturbation (introduced by substituting a B-atom 

o.~t 

for an A-atom) being localized in each atom. Consequently, 
it is imperative that an extension of the SS-CPA to include 
non-localized perturbations should be formulated to discuss 
the electronic properties of disordered binary alloys. One 
such extension is the subject of the present paper. 

We consider a tight-binding Hamiltonian 3C of a 
random binary alloy in which the atomic potentials as well 
as hopping integrals are assumed to be dependent on the 
composition of the alloy*. In the spirit of the coherent 
potential theory, an effective Hamiltonian K q is introduced 
via the diagonal coherent potential £ 0 and the off-diagonal 
coherent potential £ . The latter quantity is, however, 
restricted to a pair of nearest neighbor sites. This is an 
important assumption and essential to keeping the formalism 
tractable and the numerical work manageable. With this 
assumption, the diagonal and the off-diagonal randomness (ODR) 
can be treated separately. The presence of ODR necessarily 
correlates a given site l with any of its Z nearest neighbors. 
Consequently, such a correlation should be included in calcu- 
lating the effects of diagonal randomness. This requires 

Some aspects of the present work have been previously re- 
ported by Moorjani et al . (1971) and Tanaka et al. (1972). 



that the diagonal randomness should be treated in the pair 
approximation within the framework of the coherent potential 
theory. 

The pair scattering, in the absence of ODR, has been 
discussed in the literature by various authors (Aiyer et 
al . 1969, Freed and Cohen 1971, Cyrot-Lackmann and Ducastelle 
1971, Nickel and Krumhansl 1971, Schwartz and Siggia 1972, 
Leath 1972, Cyrot-Lackmann and Cyrot 1972, Schwartz and 
Ehrenreich 1972). Leath (1972) has pointed out the differ- 
ences which exist amongst various results. The main reason 
for these differences is discussed in Sections II and III 
which contain the general formulation of the problem. 

The detailed calculations are carried out in the two 
appendices. Section IV contains the main results of “this 
work and its relationship to previous approximations. The 
numerical results are discussed and compared with earlier 
calculations in Section V, and the conclusions are presented 
in the final section. 



- 6 - 


II. FORMALISM 

We consider the tight-binding model of a random 
binary alloy A B- which is characterized by the Hamil- 

X J.“X 

tonian, 

K = I + I U>W £m <m| (1) - 

i Jl/m 

The summation in the first term extends over all atomic 
sites while that in the second is restricted to the nearest 
neighbors only. The atomic energies e as well as the 
overlap integrals W^ m are taken to be composition depend- 
ent; assumes the values or fg depending on 

whether the site is occupied by an A-atom or a B-atom 

and W„ takes the values W. . , W or> or W A „ (=W r , A ). 

Jim AA’ BB AB BA 

The Green's function corresponding to the above” 

\ 

equation is defined by the relation, 

G(E) = — (2) . 

E-3C ’ 

where E is the complex energy. 

In the general spirit of the coherent potential 
theory, we introduce an effective medium for the motion 
of an electron, and assume that the effective Hamiltonian 

can be written as, 

\ 

: ' K 0 = I U>E<i| + \ | je>S <m| (3), 

l Vm 

where the summation conventions are the same as in Eq.(l). 

\ 

The effective or coherent potentials L (=S..) and 

0 Xj JL 
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£ (=£„ ) are as yet unknown, to be determined from an 

i v £m ■ 

appropriate self-consistent condition. We assume that the 
largest contribution to the off-diagonal coherent potential 
comes from the nearest neighbor sites; an assumption which 
is essential for keeping the numerical work within reason- 
able bounds. 

The static properties of the system are determined 
from the configurational average of the Green's function 
[Eq.(2)] over all possible configurations of the random 
alloy. The electron density of states, for example, is 
given by the well-known relationship, 

p (E) = - i Im Trace <G(E)> (4), 

where the angular brackets denote the configurational 
average. When the configurationally averaged medium is 
taken to be the effective medium defined by Eq.(3), then 
one obtains the identity, 

<G(E)> - G (E) = — L- (5). 

o • h— jv 

o 

The exact G(E) [Eq.(2)] is then related to G (E) by 

o 

the relationship 

G = G + G TG (6), 

0 0 0 ’ 

where the T-matrix is defined by, 

T = T(1 + G T) 

O 


( 7 ), 
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with r = 3C - K (8), 

o 

If one now takes the configuration average of Eq.(6), one 
obtains 


<T> = 0 


(9), 


a condition which can be used to determine the coherent 
potentials E and E . It should be noted that Eq. (9) 

o i ' 

represents a general exact condition. The various approxi- 
mations become clearer if the physical contents of the 
above mathematical formulation are made a bit more trans- 
parent . . • 


The actual potentials and W^ m which an 

electron experiences during its motion in a given disordered 
alloy are replaced in the effective medium by the unknown 
quantities E q and E^ . The multiple scattering of the 
electron are described by the T-matrix (Eq.7) and the unknown 
potentials determined from the condition that there be no 


further scattering in the effective medium. Since there are 
only two unknowns, we need just two equations; these are 
obtained by taking the diagonal and the off-diagonal matrix 
elements between nearest neighbors of Eq.(9). One thus 
obtains 


<T> 


ZJL 


0 


(10-a) 
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<T> im -0 (10-b 

At this point, a digression is essential to point out the 
relationship of the above formalism with the SS-CPA, and 
more important to point out the subtle but significant 
difference between our approach and that of Foo, et al. (1971) 
In SS-CPA one concentrates on a single effective atom 

and requires that there be no further scattering from it. 

The exact self-consistent condition (Eq.9) is thus replaced 
by the one which requires that the configurational average 
of the atomic t-matrix must vanish; this single equation 
being adequate to determine the single unknown . Foo, et 
al. (1971)carry out a straightforward extension of the 
SS-CPA, replacing a single atom by a pair of nearest neigh- 
bor atoms requiring that there be no further scattering 
from the two-atom cluster in the effective medium. This 
clearly does not treat the Z nearest neighbors of a 
given atom on the same footing; a drawback realized but not 
accounted for by Foo, et al. (1971). It is not correct, as they 
state , that such a difficulty is an essential element of 
any self-consistent treatment of pair-scattering. As 
discussed below, the difficulty can indeed be removed and 
has the expected effect of multiplying all scattering matrix 
elements from a pair of nearest neighbors by a factor of Z. 

To illustrate this point, let us decompose the T-matrix as, 
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T = y T(j&) + y T (1 ,m) + y T.(£,m,n) + (11)., 

I l/m n/(£,m) 

where, T ( i) is an operator which accounts for scattering 
"t h. 

■from the j i site, T(£,m) from all pairs of sites, 

T(£,m,n) from all three-site clusters and so on. The first 
summation is over all sites, the second over all pairs and 
so on. In., the approximation, where the effects of three 
atoms and higher order clusters are neglected and the two- 
site clusters are restricted only to the nearest neighbors, 

Eq. (11) can be truncated and written as, 


T 


£ T(jO + T(£,m) 

i m^ & 


( 12 ), 


where the second summation is only over the nearest neigh- 
bors. Taking the configurational average, we obtain 

<T> = } <T U)> + y < T ( X , m) > (13). 

Z— < L-- 

l m ^ Z 


The matrix elements in Wannier representation are therefore 
given by, 




<T bi " * L 


and 


<T> ’ = / <T(Jt,m)> . 

£m i—i ’ ' Zm 

m^£ 


(14-a) 


(14-b) 



ns ns 

Combining Eq (14) with Eq (10), and taking advantage of 
the fact that configurationally averaged quantities are 
translationally invariant, one finally obtains, 

<T(4)> M + Z<TU,m)> ^ =0 

and 

Z<TU,m)> Am =0 

Foo, et al. (1971) leave out the factor Z in the above 
equations which as seen is important only in Eq.(15-a). 

In relation to the work of various authors on pair 

scattering, this point is further discussed in the next 
section where we obtain explicit expressions for the 
diagonal and the off-diagonal matrix elements of the 


(15-a) 


(15-b ) 


T-matrix 
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III. MATRIX ELEMENTS OF THE T-MATRIX 

The form of the actual and the effective Hamiltonians 
[Eq ns (l) and (3)] suggest that the perturbation Hamiltonian 
r (Eq.8) can be conveniently separated into the diagonal 
and the off-diagonal parts, as 

r = r (1) + r (2) 

where, 

r ( D = y r ( 1 ) (i) = Y UXe - s Xi\ 

1—1 1—1 Xi 0 

J l i 

and 

r <2) * I r ( 2 ) (a) - I U> (w im - s Xm| 

(O') m/jl 

In the last equation, a. denotes a pair of nearest neigh- 
bor sites ( j£,m) . 

Corresponding to and , one can now 

define the T-matrices via the equation 

T (i) = r (l) [l + G Q T (l) j (i = 1,2) (18). 

T^^ and T^ 2 ^ describe the multiple scattering of an 
electron from the diagonal and the off-diagonal perturba- 
tion Hamiltonians respectively. The total T-matrix (Eq.7) 


(16) 


(17-a) 


(17-b) 
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can then be written in terms of and T^ 2 ^ as, 


-[ 


R 


Q + |T< 2 > + S P 


( 19 ) 


where , 


P = 


1 - G S 
L o 


Q 


= [i - c 

L o 


1 
R i 


-1 


-1 


R = 


t (2) g t (1) 


(20-a ) , 
(20-b ) , 
(20-c ) , 


and , 


S = T^G T ^ 
o 


(20-d ) . 


It should be noticed that Eq. (19) represents an exact 
expression for T in the nearest neighbor pair approxi- 
mation and includes the sum of the scatterings from F^^ 
( 2 ) 

and r plus all the correlations where an electron is 
alternatively scattered any number of times by and 


,( 2 ) 


To determine the matrix elements of T, one needs 


to calculate the matrix elements of T^^ and T^ 2 ^ only 
since P, Q, R, and S are expressed in terms of these 
quantities (Eq.20). 

(2) 

To evaluate the matrix elements of T , we com- 
bine Eq ns (17-b) and (18) to write, 


, (2 ) = V r< 2 > 


, r^'(co + ) r (2) (a)G r^ ; (|8) + .. 

t—j is 0 

a. a, (3 


i(2) 


( 21 ), 


-14- 


which can be rewritten in terms of a pair-site t-matrix, 


t (a) = r (2) (a)f 1 + g t (a) 


as , 


v? - v 1 

T ( J l t(a) + ^ t(a)G o t(/3) + 


a 


a^/3 


Thus in the 1 nearest-neighbor pair approximation, 


<t<2)> m - Z 


and . 


<T (2) > ta - Z <t im > 


where , 


t M = <4|t(a)U> and t ^ = <jj|t(a)|ra> 


The matrix elements t . . and t „ are obtained from 

SUL lm 

Eq. (22) and solving a set of coupled equations to yield, 


t JU T ml g o / i (1 " r m£ g i ) r m^ g 0 J 


and, 


t Im r to (1 “ r jem g i (1 _ r mX g i )S " T mZ g o J 


In Eq nS (27) and (28), T = T . = W „ - E and the 

’ mJL Urn Jim i 

diagonal and the off-diagonal matrix elements of G 


are given by, 


( 22 ), 


(23). 


(24), 


(25), 


(26). 


(27) 


( 28 ). 
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g = <x|g | i) and g = < | G | m) (29). 

0.0 1 o 

The presence of off-diagonal randomness , as seen in 
Eq nS (24-28), correlates the 4 th site to its nearest 
neighbor sites via r . This further manifests itself 
in the appearance of g , which determines the propagation 
of a single particle state from a given site to any of its 
Z nearest neighbors. To be consistent this correlation 
should also be included in the diagonal randomness. This 
requires that T^^ be treated in the pair approximation, 
in contrast to the recent work of Blackman, et al. (1971). 

The pair scattering treatment of T^^ has been 
carried out by various authors in the literature. How- 
ever, as Leath (1972) hasrecently pointed out, some differ- 
ences among various results still exist. A simple 
diagrammatic procedure outlined below leads to the express- 
ion first reported by Cyrot-Lackmann and Ducastelle (1971). 

Combining Eq ns (17-a) and (18), T^^ can be 

written as, 

t (1) = .^r ( 1 ) u) + £ r (1) (^ )G o r ( 1 ) u 2 ) + — (30), 



which in turn can be written in terms of the single-site 
t-matrix, 

tu) = r (1) u)[i + G o tu) ] 


(31), 
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as, 

t ( i) = ) t(D + y tu >g t(£') 

l — • i—i 10 s 

l l la • 

1 2 

+ 7 t(A )G t(l )G t(l ) + (32) 

£-i 1 0 2 0 3 

X 7^ A 
12 3 

The summation convention in the last equation implies that 
no two successive indices can be equal to each other. For 
example, in the third term i and i t i but A can 

be equal to A . Diagrammatically , Eq. (32) can be written 
as shown in Fig.(l). If one now neglects the effects of 
three-sites and higher order clusters, then one needs to keep 
only those diagrams [Fig. (1)] which involve one and two dis- 
tinct sites. Thus, in the pair approximation, some of the 
diagrams which contribute to T^^ are those shown in Fig. 
(2). Summing this series of diagrams and taking the config- 
urational average, one obtains, 

<t ( i) > = y<tu)> + y< t (x )g t u )g tu ) 

t— i. 1 O 2 0 1 

A 1 /i 

1 2 

r -|-^ 

l-Gt(A)Gt(A) + t(A )G t(A ) 

L 0 20 lJ 10 2 • 

r t! 

1 - G t(A )G t(A ) ! > (33) 

• L • 0 10 2 J 

Next, we restrict A to be the nearest neighbors of A , 

2 ' - 1 

and take diagonal and nearest neighbor off-diagonal matrix 
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elements in Wannier representation. Using the isotropy 
and the translational invariance of the lattice we find 


(T (1) > J - 


diag 


V <T (1) ) 

U W ' ZZ 


(34-a ) 


where. 


<T<1)> « " <V + Z 


(34-b) 


and, 


<T (1) > = Y <T (1) > 

off-diag Urn. 

(!,m) 


(35-a) 


with 


<T (1) >. •= Z <t.t g /(1-t.t g 2 )) 

j0m X m to i Z m to i 


(35-b ) 


where m is one of the nearest neighbors of Z and, 


= UltU)U> 


V S o 


! - g (€ E ) 


(36) 


The last equality follows from Eq.(31) combined with Eq . 
(17-a) . 

The expressions (34-b) and (35-b) are easily shown 
to be identical to the results obtained by Cyrot-Lackmann 
and Ducastelle (1971) but differ from the treatments of other 
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authors (Nickel and Krumhansl 1971, Leath 1972, Foo et al. 
1971) by the appearance of Z the number of nearest 
neighbors. It should also be pointed out that the recent 
criticism by Cyrot-Lackmann and Cyrot (1972) of their 
earlier formulation (Cyrot-Lackmann and Ducastelle 1971) is 
valid only if one considers atoms which al*e not nearest 
neighbors of each other. 

It remains to evaluate the matrix elements of P, Q, 
R and S [Eq.(20)3 and take the configurational averages. 

This part of the calculation is carried out in the two appen 
dices and the results are discussed in the next section. 
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IV. RESULTS 


As shown in Appendix II, we obtain two equations 
from the two self-consistent relations [Eq nS (10-a) and 
(10-b)] which can be written as the coupled equation for 
the two coherent potentials. The resulting equations, 
obtained from Eq.(A-23) and the combination of Eq. (A-29) 
with Eq. (A-20) , respectively are 


£ = U - (€.- £ )g (c — £ ) + Z 1 - (e - £ )g i 

o o A o o B o L A o o J 


t 1 - <V S o>8 0 > o (37) 


and, 


= U 1 + ^ + L 1 " Jl 1 - 


[l-.g i r(y,y>] (f i+ fJ 


(38) 


The quantities F , F , F , £ and T's are defined by Eq nS 

0 12 

(A-17), (A-27 ) , (A-21) , (A-30) and (A-19-f) respectively. 

U q and U are the potentials in the virtual crystal 
approximation given by, 


U o ■ xt A + yE B <39> 

and 

- * S W AA + 2xyff AB + ySff BB 


(40) 
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One notices that previously used . approximations 
ns 

are contained in Eq (37) and (38). Thus the first terms 
in these equations represent the virtual crystal approxi- 
mation and the first two terms in Eq.(37) reproduce the 
result obtained in the SS-CPA. The rest of the terms are 
the corrections due to the inclusion of pair-scattering and 

off-diagonal randomness in the present analysis. 

ns 

To solve Eq (37) and (38), it is convenient to 
transform them in dimensionless form. We shall therefore 
refer all quantities to ZW the half-bandwidth of the 


(41) 


(42) 

(43) 

(44) 

\ . . 

\ • ’ • • 


nil 

pure A-lattice. Furthermore, we assume that e = - €., 

A 

and obtain 


o - = (x - y) 6 + (6 2 - a 3 ) p 


O ZW 


AA 


o o 


ir 


+ Zj 1 - (6 - a )P || 1 + (6 + a )p if 

L 00-IL OOJ 


and, 


where , 


cr = — — = W + £'+ (1 - T. P ) (1 - r._ P ) 
l zw AA 1 AA 1 AB 1 


V 


V- V ZW AA' «'-« /ZW AA 


P = ZW..g 
o AA o 
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P x - ZW AA^i 


(45) 


5 - ^A" €b - c /™ 
o c a /zw aa 


2ZW 


AA 


(46) 


W = — (x 2 + 2yxy + 0y 2 ) 


(47) 


and , 


r., - r(x ’ x) -i-o 

AA ZW AA 2 1 


t ~ t _ ^( x > y ) 
AB zw AA 


- a 


r = r ^y.y) = £ 

BB 


ZW 


AA 


(48) 


The last equations result from defining the bandwidths of 

A-B and B-B lattice in terms of the A-A lattice; 

W . = yW AA and A . . 

AB AA BB AA 

Before representing the numerical results, let us 

note that g and g (Eq.28) can be written in k-repre- 
o 1 

sentation as, 


= 1 V 1 

0 N f E - S - L y(k) 
k o i 


(49) 


and 


e , — 


l V 


ik . (I -m) 


u 


N ~ E - L - E y (k) 

K ox 


(50) 
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where y(k) = > is the structure factor, £ being 

L i 

t 

the nearest-neighbor vector. In Eq. (50) , i and m are 
nearest neighbors and therefore, 



(51) 


where the dimensionless parameter ot is defined by, 


E - £ E' - CT 


a. = 


ZE 


Za 


(E ' = E/ZW aa ) 


(52) 


The dimensionless quantities P° and P^ [Eq (44) and 
(45)] are therefore given by 


P = — 
o ZCT, 


1 V 

N L 


a - 


z - 1 y (k) 


(53), 


and 


(y. P — 

o ZCT. 


(54). 


The quantity of interest is the electronic density of states 
(Eq.4) which can be expressed as, 


P( E) - - i Im E <l|G U> = - - Img 

77 . O 77 


or 


P ' (E) = ZW 


AA 


(Pi J 
' N 


(E)Nl = - - ImP 
77 0 


(55), 
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p'(E) being the density of states/a.tom expressed in units 

of half-bandwidth of pure A-lattice. 

ns 

The two non-linear Eq (41) and (42) are now 

solved iteratively for a and a as a function of 

.0 1 

reduced energy E' , for given values of x, the concen- 
tration of A-atoras; 6, the potential difference between 
type-A and type-B atoms, and /3 and y the relative band- 
width of B-B and A-B lattices respectively with respect 
to the A-A lattice. The knowledge of the coherent poten- 
tials and a allows one to compute the imaginary 

part of P and the density of states is obtained from 

o 

Eq. (55) . We have carried out the numerical calculations 
for a number of different alloys for the body-centered 
cubic lattice and the results are discussed in the next 


section . 
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DISCUSSION 

The numerical results for electron density of 
states/atom [in fact n p', where p' is given by Eq.(55)] 
as a function of energy for various alloy parameters are 
presented in this . section. The results are illustrated by 
a number of figures where the normalized Energy E' (Eq. 
52) has been shifted by 6 so that the band corresponding 
to the pure A-lattice in these dimensionless units is 
located at the origin and extends from - 1 to +1. 
Throughout the present analysis we have considered the 
body-centered cubic (BCC) lattice for which the density of 
states was obtained by using the lattice Green's function 
recently evaluated by Morita and Horiguchi (1971). The 
density of states for ordered pure A-BCC lattice is shown 
in Figure 3 and all effects of disorder are measured with 
respect to it. Since we have chosen = “ the 

effect of substituting B-atoms for A-atoms is expected to 
be the appearance of a "minority band" located at -26, 
whose width is controlled by the parameters x, /?, and y. 
Whether any structure due to presence of minority atoms is 
discernible in the density of states curves, depends as 
shown below, on the values of these parameters and 6. 

j The results for the case of a dilute alloy 
(x = 0.9) are presented in Figures 4 and 5 for 6 = 0.2 

\ • ' ■ 
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and two different values for the set (/ S , y ) . As shown in 
the figures, the general effect of disorder is to remove 
the singularity at the origin in the density of states of 
the ordered lattice and modify the bandwidth, which is 
increased for the alloy represented by (yS , y) = (1.4, 1.2) 
[pig. 4] and decreased for the case (£, y) = (0.15, 0.1) 
[pig. 5]. The density of states in the vicinity of 
Energy = - 0.4 is increased relative to the ordered 
lattice showing the presence of the contribution of the 
minority band. For the larger bandwidth case [Fig. 4] the 
minority band is relatively flat and no visible structure 
shows up at Energy = - 0.4. This effect, as expected, is 
more pronounced for the shorter bandwidth alloy and, as 
shown in the insert of Figure 5, a shoulder in the density 
of states is indeed discernable. For a non-dilute alloy, 
still more pronounced structure would be expected and this 
is indicated in Figures (6)~(8) for x = 0.6, 6 = 0.2 and 
three different bandwidth combinations (£, y). Note that 
the peak due to the minority band does not lie exactly at 
-0.4 but slowly shifts towards higher energy with increas- 
ing values of fi and y, while the main peak shifts towards 
lower energy. The magnitude of the impurity band peak is 
lowest for (/8, y) = (1.4, 1.2) and highest for 
( 0 , 7) = (0.4, 0.7). 
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In the cases discussed above, we have only 
considered alloys where the main band and the minority 
band overlap in energy (6 < 1). In order to compare 
results with earlier calculations, it is desirable to 
consider the split-band case (6 > 1) so that the effects 
of the minority band are clearly separated from that of 
the main band. Also to facilitate comparison with the 
SS-CPA, we first consider the alloy with no off-diagonal 
randomness, i.e. /3 •= y = i. For this particular case, the 
results of Blackman, et al. (1971) reduce to the SS-CPA. This 
is due to the neglect of pair scattering in their work. 

Our results (labeled TS-CPA for two sites-CPA) for x = 0.9 
and 6 = 1.05 are presented in Figure 9 and compared with 
the SS-CPA.. In the main band, the differences between 

the TS-CPA and SS-CPA results are negligible and for clar- 
ity only TS-CPA results are shown. However, important and 
major differences appear in the shape of the minority band 
(Fig. 9) as obtained in the present work and that calcu- 
lated from SS-CPA. In the later, one obtains a smoothly 

varying curve for the density of states in the minority 
band; a result which originates from the fact that in the 
SS-CPA, the perturbation is localized at a given site. 
Consequently the coherent potential is k-independent and 

I 

is incapable of sampling the lattice. The k-dependent 
coherent potentials in the present formulation lead to the 

\ ' * ; 



/ 
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appearance of structure on both sides of the minority band 
(Fig. 9). Following quite a different formulation, Schwartz 
and Siggia (1972) alsofind structure in the minority band. 

Foo, et al. (1971) in their calculations for one-dimensional 
random chain, have also reported structure in the minority 
band; however, the location of their minority band is a 
strong function of the bandwidths /S and y. This incorrect 
result is due to their special choices of the values of £ 
and y and is typical of the shape of the density of states 
of one-dimensional chains. 

It is of interest to plot the real and imaginary 
parts of the coherent potentials and cr^ and compare 

them with the SS-CPA. This is done in Figures 10-12. 

As shown in Figures 10 and 11, the inclusion of pair 
scattering leads to the presence of structure in ReCT o and 
Ima which is responsible for the structure in the density 
of states (Fig. 9). For the off-diagonal coherent potential 
(Fig. 12), the SS-CPA results would be Imcr^ = 0 and 
Reo^ = 0.125 over the entire range of the minority band. 

Finally, Figure 13 shows the density of states for 

the case when off-diagonal randomness is present. The 

c 

presence of the off-diagonal randomness has the advantage 
of controlling the bandwidth of the minority band by 
varying the values of /S and y. As shown in Figure 13 
for the case y = 0.5, further structure appears in the 
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density of states of the minority band due to the splitting 
of the central peak. The structure in the absence of the 
off-diagonal randomness (Fig. 9) has been assigned by 
Schwartz and Siggia (1972) to the bonding and anti-bonding 
levels of a molecule embedded in an effective medium. At 
the time of the present writing, we have not found a satis- 
factory explanation for the extra structure seen in the 
presence of the off-diagonal randomness (Fig. 13). 
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SUMMARY 

• In this paper, we have extended the single-site 
coherent potential approximation for random binary alloys 
to include the effects of the off-diagonal randomness and 
the pair scattering. The formalism reduces to previous 
approximations used in this problem in appropriate limits. 
The numerical results for a number of different alloys are 
reported and it is shown that the inclusion of pair scatter- 
ing and off-diagonal randomness leads to the appearance of 
structure in the density of states of the minority compon- 
ent band. 

There remains the problem of investigating the 
effect of scattering from clusters higher than the two-sites 
clusters. In principle, this is quite straightforward in 
the framework of the present theory; however, the numerical 
work involved is expected to be quite formidable. The 
.theory can also be applied to the calculation of conduct- 
ivity (Sokolo'ski et al. 1972) and indeed the study of 
disordered magnetic systems (Bose et al. 1972). The 
detailed calculations of such extensions will be reported 


in the near future. 
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APPENDIX I 


The Matrix Elements 


The matrix elements of R [Eq.(20-c)] and S [Eq. 
(20-d)] between Wannier states on a nearest neighbor pair 
of sites (i,m) are easily written down in terms of the 


matrix elements of T 


( 1 ) 


and T 


( 2 ) 


as 


R ii R im 


R . R 
mi. mm. 


.(2) rp (2)' 

ii 


im 


( 2 ) ( 2 ) 

mi mm 


g 0 g ! 


g i g o 



r (1) 

1 im 

r (D 

mm 


(A-l) , 


and , 


S ii S im 


S „ S 
mi mm 


T d) T (l) 

ii A im 

T d) T (l) 

mi 1 mm 


g o g x 


g x g o 



,( 2 ) „ ( 2 ) 


ii 


mi 


im 


, ( 2 ) ^,( 2 ) 


mm 


CA-2) 


We notice that, 

S . . = R . . , S. = R . , S . = R . , and S = R (A-3) 

ii ii’ im mi’ mi im’ mm mm 

The matrix elements of P and Q can now be obtained by 
taking the diagonal and the off-diagonal matrix elements 
of Eq. (20-a) and Eq. (20-b)', the resulting expressions are, • 


Q ii Q im 


Q „ Q 

m i \m , 


1 - g R_ - g R. g R. + g R m 
o mm i im o im i mm 


D(i,m) 


g o R mi + Mm 


1 g o R i< _ g i E n<. 


and , 


(A-4 ) 
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P P 
ii lm 


P , P 
mi mm , 


D(4,m) 


^ - ^o^iran g i^im g o^£m + ^l^mra 
g o S m« + * t S ll 1 - g o S U 6 , S .l ' (A - 5) ’ 


where , 


D(£,m) 1 - S 0 ( R ££ + R mm ^ “ g i ^ R jem + R m^ + 


(g o - O (R M R mm - R WW 


(A-6) 


ns 

The same D(4,m) appears in Eq (A-4) and (A-5) due to 
the symmetry relations between R's and S's [Eq.(A-3)]. 

It is now straightforward to write down the matrix 
elements of T (Eq.19); we, however, follow a slightly 
modified approach in order to get the final equations in 
the form which can be easily compared with previous calcu- 
lations. To evaluate the diagonal matrix elements T we 

separate out the single-site and two-site contributions, 
leading to 


T = t . + Z T , „ ( & ,m) 
li l il 


(A-7 ) 


where , 


T „(J!, ra ) 


<T « )+ + (T ta + R l J\l 


+ (T il >+ S H )P M + (T im >+ S lm )P n,l 


(A-8), 


with 
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T W > - - VmO' 


(A-9 ) 


T lm ~ t « t .n E i /<1 - VmO- 


rp( 2 ) = t and T^ 2 ^ = t 

zA’ T fm 1m 


(A-10) 

(A-ll) 


The quantities t.. and t. in Eq. (A-ll) are given by 

Jo Xj -oHl 

Eq. (27) and Eq. (28) respectively. . 

To evaluate the off-diagonal matrix elements, we 

write 


- Z 


£m 


( A-12 ) 


where 


W*- m) 


T -- )+ ( T M 


im 


U + «£*>>«*» + ( T im >+ mm 


, T T 1 4- £T T ^ 2 ^ 4- V T ^ ^ 1 S fa T ^ 2 ^ 4 cr T ^ ^ I C! Ip 

iL 1 + S o T JU + S i T Jhn J S ff ii + g o T M J S m£j P £m 


fr (2) (2)1 

+ u 1 + g o T L + g , T L J s , m + 


I g t< 2) + g IS Tp 

L°i 1 1 1 Zm J mmJ mm 


(A-13 ) 


The last equation is obtained by rewriting Eq.(19) as, 


_ rp( 2 ) ,„,(!) 


,( 2 ), 


T = T '+ (T v y + R)Q + (S + T v 'G S)P (A-14) 


It remains to take configurational average of T.. and 

AJ fj 


T^ m and this is carried out in the next appendix, 
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APPENDIX II 

Configurational Averages 

Taking the configurational average of Eq. (A-7) and 
invoking the self-consistency condition [Eq.(lO-a)], we 
obtain 


<t A >-'.= xt x + yt y = - Z<T AA (/,m)> . 


where, from Eq. (36) , 


f A,B “ S o 


x,y i - g (e - S ) 
& o K A, B o' 


Also , 


<T jeA ( f,m)) 


V" 

F = / 

o Z_, 


IJ t m ( i ’ j) 


I=x,y 

J=x,y 


with, 






r ( -I \ 1 

+ L T k (1 ' J) + 


+ + s «0,J)]p M (I,J) 




(A-15), 


(A-16 ) . 


( A-17 ) 


(A-18 ) 


\ 



-35- 


R (I,J) etc. are obtained from Eq . (A-l) , (A-2), (A-4) 

AJ AJ 

and (A-5) and various T^^'s and T^^'s from Eq 11S 
(A-9)- (A-ll) . One then has, 


t 2 t 2 2 

TW(I.J) - — 

1 - 


(A-19a ) , 


I J S i 


(A-19b) , 


Ttt'lM) - t3jt;g ' , 

ram 1 - vx 


(A-19c ) , 


(2) 


_ a,(2) 


T ]/ (I ’ J) = mm ' ( 1 > J ) 


r 2 d,j)g r 


[i - r<i, j) gi ] 2 - r 3 (i,j) g = 


(A-19d) , 


and 






r(i, j)[i - r(i, j) gi ] 

Ci-r(i,j) gi ] 2 - r 2 (i,j ) g 2 


(A-19e ) 


where 


T(i,j) = W(I,J) - 2^, ( W(x , x ) = W AA , etc.) (A-19f ) 

Finally taking the configurational average of Eq. (A-12) 
and invoking the self-consistency condition [Eq.(lO-b)], 
we obtain 

< T } > = " <F > 

Zxn 2 


(A-20) 
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. where, 

<F > = Y IJ F (I,J) (A-21 ) 

2 f-< 2 

I-x,y 

J=x,y 

and F (I.J) is obtained from Eq.(A-13), 

2 

+ {L 1 +So T « )(I ' J)+g 1 T im )( I- J, ] S « (I ’ J) 

+ [ g i T ii )<I , J) +g o T im ) < I - J )] S m^ (I - J )} P im a - J) 

+ IC 1 + g o T ii >(I ’ J) + g 1 T im )(I ’ J) ] S tm <I ’ J) 


(A-22) 


The two equations of interest are Eq. (A-15) and 
Eq.(A-20) and to put them in reasonable form, we first 
rewrite Eq. (A-15) as, 


€ . — E 

A o 


1 - g o<V S o> 


+ y 


€ - E 

B o 


= - Z<T ( f ,ra) > = - ZF 


ll 


1 - g o <C B‘ S 0 5 
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or, 


£ = U - (c - £ )g (e - £ ) + Z 1 - (e - £ )g j 

o o A ooB o L A ooJ 




1 - (€ — £ )g jF 

fc> O O J 0 


(A-23 ) , 


which is Eq.(37). 

To obtain Eq.(38), we rewrite 


-.( 2 ) = 

' im 


im 


i - r. g 


+ F 


(A-24), 


where , 


F = 
i 


•p 3 2 

r Am g o 


(i - r. g )Ci - r. g ) 2 - rf g 2 ] 

v £m fc> i 7 jhn* 3 ! £m & o 


(A-25) . 


Therefore, 

<T ( , 2) > - x S r(x ’ x) + 2xy r . < X . ; . y i— 

m l-g 1 T(x,x) l-g^TCx.y) 

+ y 2 — ■ IlLill + <F > (A-26) , 

l-g^Cy.y) 1 

where, 

(F^ = ^ I J F ( I , J) (A-27 ) , 

I=x,y 
J-x , y 

i . 

i 

and , 

\ 

\ * ' ■' 
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r 3 u, j)g 2 

F (I,J) = 2 (A-28) . 

[i -r(i,j) gi ][ (i -r(i J j)g i ) 2 -r s (i,j) g ^] 

From Eq.(A-26), we have, 

U - £ + i 

_ i i 

Cl-g i r(x,x)][l-g i r(x,y)][l-g i r(y,y)] 

+ <F > (A-29 ) , 

the value of the overlap integral in the virtual 
crystal approximation is given by Eq.(40), and 

i = r(x,x)r(x,y)r(y,y)gj- (x 2 + 2xy) 

r(x,x)T(x,y)g i - (x 2 + y 2 )r(x,x)r(y,y) gi 

- (y 2 + 2xy )F(x , y )T(y , y ) g^ (A-30) . 

Combining Eq. (A-29) with Eq. (A-20), we finally obtain 
Eq. (38) for 2^ . 


<t (2) > 


where U , 
i 
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FIGURE CAPTIONS 

Figure 1. The diagrammatic representation of the first four 
terms in Eq. (23) . The dashed vertical line repre 
sents scattering from a given site denoted by a 

star. The horizontal line (G ) represents the 

o 

motion in the effective medium. All sites are- 
distinct and summation over distinct sites is 
implied . 

Figure 2. The first diagram represents the scattering from 
single-site clusters while the others are some 
of the diagrams obtained from Fig. (1) which are . 
appropriate to scattering from two-site clusters. 

Figure 3. - The density of states as a function of energy for 

an ordered body-centered cubic lattice. 

Figure 4. The density of. states Vs. energy for a random 
A 0.9 B 0.1 alloy; > 1 ’ 

Figure 5. The density of states Vs. energy for a random 
A 0 9 B 0 1 a H°y; /S , y < 1. The insert shows 
the density of states in the vicinity of the 
energy where one would expect the effects aris- 
ing from the presence of B-atoms. 

Figure 6. The density of states Vs. energy for a random 
A 0.6 B 0.4 


alloy; /3 , y > 1. 
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Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


7. The density of states Vs. energy for a random 

A 0.6 B 0.4 alloy; y < 1 * 

8. Same as for Fig. 7 but diffez'ent values of /8 
and y. 

9. The density of states Vs. energy for a random 
A„ n B n .. alloy in the split-band case (6 > 1) . 
Note the break in the horizontal scale. The 
vertical scale on the right is for the main band 
while that on the left is for the "minority!' 
band. The results of the present work are 
labeled TS-CPA. 

10. The real part of the diagonal coherent potential 

a over the energy range of the "minority" 
o 

band. Comparison between present results 
(TS-CPA) and SS-CPA. 

11. The imaginary part of the diagonal coherent 
potential a q over the energy range of the 
"minority" band. Comparison between present 
results (TS-CPA) and SS-CPA. Note the break 
in the vertical scale. 

12. The real and the imaginary part of the off- 
diagonal coherent potential a over the 
energy range of the "minority" band. 



Figure 13. . The density of states Vs. Energy for a random 

Aq g B 0 ^ alloy in the split-band case (6 > 1) , 
including the off-diagonal randomness. Note 
the break in the horizontal scale. The vertical 
scale on the right is for the main band while 
that on the left is for the ’’minority” band. 
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